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Abstract. It is hard to imagine curved spacetimes of General Relativity. A simple but powerful way how to achieve this 
is visualizing them via embedding diagrams of both ordinary geometry and optical reference geometry. They facilitate 
to gain an intuitive insight into the gravitational field rendered into a curved spacetime, and to assess the influence of 
parameters like electric charge and spin of a black hole, magnetic field or cosmological constant. Optical reference geometry 
and related inertial forces and their relationship to embedding diagrams are particularly useful for investigation of test 
particles motion. Embedding diagrams of static and spherically symmetric, or stationary and axially symmetric black-hole 
and naked-singularity spacetimes thus present a useful concept for intuitive understanding of these spacetimes' nature. We 
concentrate on general way of embedding into 3-dimensional Euclidean space, and give a set of illustrative examples. 
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INTRODUCTION 

The analysis of embedding diagrams GlQ>liLliLS>S>l2lS] rank among the most fundamental techniques that enable 
understanding phenomena present in extremely strong gravitational fields of black holes and other compact objects. 
The influence of charge, spin, magnetic field or cosmological constant on the structure of spacetimes can suitably be 
demonstrated by embedding diagrams of 2-dimensional sections of the ordinary geometry (t — const hypersurfaces) 
into 3-dimensional Euclidean geometry and their relation to the well known 'Schwarzschild throat' Q9|] . 

Properties of the motion of both massive and massless test particles can be properly understood in the framework of 
optical reference geometry allowing introduction of the concept of gravitational and inertial forces in the framework of 
general relativity in a natural way lUOl [ill Il2l U3l \\A. U5l UOL U7L U8I1 . providing a description of relativistic dynamics 
in accord with Newtonian intuition. 

The optical geometry results from an appropriate conformal (3+1) splitting, reflecting certain hidden properties of 
the spacetimes under consideration through their geodesic structure. Let us recall the basic properties of the optical 
geometry. The geodesies of the optical geometry related to static spacetimes coincide with trajectories of light, thus 
being 'optically straight' fnl [l^l . Moreover, the geodesies are 'dynamically straight,' because test particles moving 
along them are held by a velocity-independent force 111 211 ; they are also 'inertially straight,' because gyroscopes carried 
along them do not precess along the direction of motion lfl3ll . 

Some fundamental properties of the optical geometry can be appropriately demonstrated namely by embedding 
diagrams of its representative sections flol I3I Ejj. Because we are familiar with the Euclidean space, 2-dimensional 
sections of the optical space are usually embedded into the 3-dimensional Euclidean space. In spherically symmetric 
static spacetimes, the central planes are the most convenient for embedding - with no loss of generality one can 
consider the equatorial plane, choosing the coordinate system such that — it /2. In Kerr-Newman backgrounds, 
the most representative section is the equatorial plane, which is the symmetry plane. This plane is also of great 
astrophysical importance, particularly in connection with the theory of accretion disks @L 

In spherically symmetric spacetimes (Schwarzschild [17], Schwarzschild-de Sitter [1], Reissner-Nordstrom 01, 
and Reissner-Nordstrom-de Sitter [4]), an interesting coincidence appears: the turning points of the central-plane 
embedding diagrams of the optical space and the photon circular orbits are located exactly at the radii where the 
centrifugal force, related to the optical space, vanishes and reverses sign. The same conclusion holds for the Ernst 
spacetime ]5[]. 

However, in the rotating black-hole backgrounds, the centrifugal force does not vanish at the radii of photon circular 
orbits f2(il, l2lll . Of course, the same statement is true if these rotating backgrounds carry a non-zero electric charge. 
Throughout the paper, geometrical units (c = G — 1) are used. 



OPTICAL GEOMETRY AND INERTIAL FORCES 



The notions of the optical reference geometry and the related inertial forces are convenient for spacetimes with 
symmetries, particularly for stationary (static) and axisymmetric (spherically symmetric) ones. However, they can 
be introduced for a general spacetime lacking any symmetry II 1 811 - 

Assuming a hypersurface globally orthogonal to a timelike unit vector field n and a scalar field <t> satisfying the 
conditions n K Vxn M — 0, n K n K — —1, hi = n K V K n / t = V^cD, the 4-velocity u of a test particle of rest mass m can be 
uniquely decomposed as 

u K = y(n K +vz K ). (1) 

Here, r is a unit vector orthogonal to n, v is the speed, and y = (1 — u 2 ) _1//2 . 

Introducing, based on the approach of Abramowicz et al. [18], a projected 3-space orthogonal to n with the positive 
definite metric h K x giving the so-called ordinary projected geometry, and the optical geometry h K \ by conformal 
rescaling 

h K x — t~ 2 ®h K x , h K x= g K x + n K nx, (2) 

the projection of the 4-acceleration — u^W^ux can be uniquely decomposed into terms proportional to the zeroth, 
first, and second powers of o, respectively, and the velocity change v = (e^y o) ijU w". Thus, we arrive at a covariant 
definition of inertial forces analogous to the Newtonian physics lfl8l. l22ll 

mat = G K (o°) + C K (o 1 ) + Z K (v 2 ) + E K (») , (3) 

where the terms 

G K = -mV K <b = -m® K , C K = -my 2 vn A (VxT K -V K Tz), Z K — -m(y v) 2 t X Vxt k , E K = -m6i K , (4) 

correspond to the gravitational, Coriolis-Lense-Thirring, centrifugal and Euler force, respectively. Here, ? is the unit 
vector along % in the optical geometry and V is the covariant derivative with respect to the optical geometry. 

In the simple case of static spacetimes with a field of timelike Killing vectors n = f if. — 8/8t, we are dealing with 
space components that we will denote by Latin indices in the following. The metric coefficients of the optical reference 
geometry are given by the formula 

r — 2<t>; 20 £>n kV zc\ 

«i<r = e h ik , e =—gnv<Z(ffi(t) = —gtt- (5) 
In the optical geometry, we can define a 3-momentum of a test particle and a 3-force acting on the particle lfl8il23ll 

/5'=e 2 V, fi=e 2 *fi, (6) 

where p l , fj are space components of the 4-momentum p and the 4-force /. Instead of the full spacetime form, the 
equation of motion takes the following form in the optical geometry, 

mf i =p k V k p i + \m 2 %0, (7) 

where represents the covariant derivative components with respect to the optical geometry. We can see directly 
that photon trajectories (m — 0) are geodesies of the optical geometry. The first term on the right-hand side of (|7]i 
corresponds to the centrifugal force, the second one corresponds to the gravitational force. 

In the case of Kerr-Newman spacetimes, detailed analysis (see, e.g., 112 ill ) leads to decomposition into gravitational, 
Coriolis-Lense-Thirring and centrifugal forces 

r ! - l"i ( S^-Sttg^ W ,„ 
G K - -m® K - -m-8 K lnl I , (8) 

2 |_ \ SH /J 

C K — mA 2 (D. - ^lnrf)VI^[ 5 k (sttfig^ 2 ^ + ^lnrf Qk-Ju 



(9) 

_I m A 2 (£2-Q LNRF ) 2 W sJlnf^^ )1 , (10) 



respectively; we also denote i — e~®n. The Euler force will appear for Q, ^ const only, being determined by Q, — 
m^V^Q (see Hill ). Here we shall concentrate on the inertial forces acting on the motion in the equatorial plane 
(9 — n /2) and their relation to the embedding diagram of the equatorial plane of the optical geometry. 

Clearly, by definition, the gravitational force is independent of the orbiting particle velocity. On the other hand, both 
the Coriolis-Lense-Thirring and the centrifugal force vanish (at any r) if Q. — Qlnrf, i-e., if the orbiting particle is 
stationary at the LNRF located at the radius of the circular orbit; this fact clearly illustrates that the LNRF are properly 
chosen for the definition of the optical geometry and inertial forces in accord with Newtonian intuition. Moreover, 
both forces vanish at some radii independently of Q. In the case of the centrifugal force, namely this property will be 
imprinted into the structure of the embedding diagrams. 



EMBEDDING DIAGRAMS 

The properties of the (optical reference) geometry can conveniently be represented by embedding of the equatorial 
(symmetry) plane into the 3-dimensional Euclidean space with line element expressed in the cylindrical coordinates 
(p , z, a) in the standard form. The embedding diagram is characterized by the embedding formula z — z (p) determin- 
ing a surface in the Euclidean space with the line element 



d^ E) = l + ( — ) \dp 2 + p 2 da 2 (11) 



isometric to the 2-dimensional equatorial plane of the ordinary or the optical space line element [2' 

de 2 = h rr dr 2 + h^d4> 2 . (12) 

The azimuthal coordinates can be identified (a = (f>) which immediately leads to p 2 — hss, and the embedding formula 
is governed by the relation 

(0) 2 =*"(S 2 - 1 - <"> 

It is convenient to transfer the embedding formula into a parametric form z(p) — z(r(p)) with r being the parameter. 
Then 

^ = ± Jhrr-i -f I • (14) 



dr 



The sign in this formula is irrelevant, leading to isometric surfaces. Because dz/dp — (dz/dr)(dr/dp), the turning 
points of the embedding diagram, giving its throats and bellies, are determined by the condition dp/dr — 0. The 
embedding diagram can be constructed if the reality condition h rr — (dp /dr) 2 > holds. 



EXAMPLES OF EMBEDDING DIAGRAMS 

In this section we present graphic representations of embedding diagrams for various spacetimes with their brief 
description. For full treatment, see the reference provided at the beginning of each subsection. 



Schwarzschild spacetime with nonzero cosmological constant 

The influence of cosmological constant on the structure of Schwarzschild-de Sitter spacetime is treated in full 
technical details by Stuchlflc and Hledflc [1]. Putting y = AM 2 /3, t —> t/M, r — > r/M, where M is the mass parameter 
and A denotes the cosmological constant, the line element in standard Schwarzschild coordinates reads 

ds 2 = -^l---yr 2 Jdt 2 +^l---yr 2 J dr 2 + r 2 (d6 2 + sin 2 6»d^ 2 ) . (15) 



Embedding diagrams of the ordinary geometry are given by the formula z — z(r), which can be obtained by 
integrating the relation 

dz _ / 2 + yr 3 \ 1/2 
dr \r— 2 — yr 3 / 

In the case of Schwarzschild-de Sitter spacetimes, the embedding can be constructed for complete static regions 
between the black-hole (r^) and cosmological (r c ) horizons. Recall that the static region exists for y < 1 /27 only. In 
the case of Schwarzschild-anti-de Sitter spacetimes, the static region extends from the black-hole horizon to infinity. 
However, we can see directly from Eq. ( fToT l that the embedding diagrams of the ordinary space can be constructed in 
a limited part of the static region, located between the black-hole horizon and r e ( rd) = (— 2/y) 1//3 . 

Embedding diagrams of optical reference geometry are given by parametric formulas z — z(r),p — p(r): 



dz 



dr r— 1 — yr 3 



4r-9-yr 4 
r(r — 2 — yr 3 ) 



1/2 / n x -l/2 

(17) 



The embedding formula z — z(p) can then be constructed by a numerical procedure. Further, it can be shown (H that 
'turning radii' of the embedding diagrams are given by the condition dp/dr = 0. Since 

we can see that the turning radius determining a throat of the embedding diagram of the optical geometry is located 
just at r — 3, corresponding to the radius of the photon circular orbit; it is exactly the same result as that obtained in 
the pure Schwarzschild case. The radius of the photon circular orbit is important from the dynamical point of view, 
because the centrifugal force related to the optical geometry reverses its sign there 1 171 12311 . Above the photon circular 
orbit, the dynamics is qualitatively Newtonian with the centrifugal force directed towards increasing r. However, at 
r = 3, the centrifugal force vanish, and at r < 3M it is directed towards decreasing r. The photon circular orbit, the 
throat of the embedding diagram of the optical geometry (dp/dr — 0), and vanishing centrifugal force, all appear at 
the radius r = 3. 

In the case of Schwarzschild-de Sitter spacetimes containing the static region, the embeddability of the optical 
geometry is restricted both from below, and from above. Using a numerical procedure, the embedding diagrams are 
constructed for the same values of y as in the case of the ordinary geometry. 

In the case of Schwarzschild-anti-de Sitter spacetimes, the embeddability of the optical geometry is restricted from 
below again; with y — > — oo the limit shifts to r — > 0, along with the radius of the black-hole horizon. The embedding 
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FIGURE 1. Left panel: Qualitative features of the embedding diagrams of the ordinary geometry of the Schwarzschild-de Sitter 
and Schwarzschild-anti-de Sitter spacetimes in a log-log diagram. One can immediately see how the diagrams with y jfe 'peel off 
the pure Schwarzschild diagram (y = 0, bold curve). All sections with y are complete (i.e., the maximum that can be embedded 
into Euclidean space is shown), except uninteresting lower parts of the throats. The diagrams clearly indicate modifications of the 
spacetime structure caused by the presence of a cosmological constant. Middle panel: Same as left panel but for optical geometry. 
Right panel: Example of embedding diagrams of the optical reference geometry of the Schwarzschild-de Sitter spacetimes. The 
pure Schwarzschild case is taken for comparison in (a). The diagrams are given for (b) y — 10 -6 , (c) v = 0.002, and (d) y — 0.03. 
They are similar to the pure Schwarzschild case, because the region near the cosmological horizon being of highly different character 
is 'cut off' by the limit of embeddability. Taken from Q]. 



diagrams are constructed by the numerical procedure for the same values of y as for the ordinary space. These diagrams 
have a special property, not present for the embedding diagrams in the other cases. Namely, they cover whole the 
asymptotic part of the Schwarzschild-anti-de Sitter spacetime, but in a restricted part of the Euclidean space. This is 
clear from the asymptotic behavior of p(r). For r — > +00, there is p ~ (— y) -1 / 2 . Clearly, with decreasing attractive 
cosmological constant the embedding diagram is deformed with increasing intensity. The circles of r — const are 
concentrated with an increasing density around p — (— y) -1 / 2 as r — > 00. 

Basic features of the embedding diagrams of the Schwarzschild-de Sitter spacetimes are illustrated in Figure[T] 



Interior uniform-density Schwarzschild spacetime with nonzero cosmological constant 

The influence of the cosmological constant on the structure of interior uniform-density Schwarzschild spacetime is 
treated in full technical details by Stuchlfk et al. 0]. 

The line element in standard Schwarzschild coordinates reads 

ds 2 = -e 2 * (r > dt 2 + e 2 ™ dr 2 + r 2 {d0 2 + sin 2 #dp 2 ) , (19) 

where 

„*« = ( X _ !f) " 1/2 , e «W = Ae -*« _ fle -* W 1 = ™ + A A = 9M = 3M-AR* 

V a 2 ) a 2 R 3 3 6M + A/? 3 6M + A/? 3 

Here M denotes total mass, p (uniform) mass density and R external radius of the configuration; A is the cosmological 
constant. Examples of embedding diagrams (for repulsive cosmological constant) are shown in Figure|2] 




Ordinary geometry Optical geometry 



FIGURE 2. Axial sections of the Schwarzschild-de Sitter (A > 0) case. The dashed curve represents pure Schwarzschild 
spacetime, the dots mark the position of the edge of the internal spacetime (i.e., r — R = xM). The embedding diagrams of 
the optical reference geometry show the existence of internal circular null geodesies which play a key role for the possibility of 
neutrino trapping. 



Reissner-Nordstrom spacetime with nonzero cosmological constant 

This section is based on Il25ll . where full-detail treatment can be found. Putting y = iAM 2 , e = Q/M, t — > t/M, 
r — > r/M, where M is the mass parameter, Q is the electric charge and A is the cosmological constant, the line element 
in standard Schwarzschild coordinates reads 

ds 2 = -^l-^ + ^-yr 2 ^dr 2 +^l-^ + ^-yr 2 ^ dr 2 + r 2 (d<9 2 + sin 2 <9 d</> 2 ) . (21) 

An example of classification based on properties of embedding diagrams of optical reference geometry is in Figure|3] 




Ernst spacetime 

The static Ernst spacetime I26I I5I1 is the only exact solution of Einstein's equations known to represent the spacetime 
of a spherically symmetric massive body or black hole of mass M immersed in an otherwise homogeneous magnetic 
field. If the magnetic field disappears, the geometry simplifies to the Schwarzschild geometry. Therefore, sometimes 
the Ernst spacetime is called magnetized Schwarzschild spacetime. Usually it is believed that for extended structures 
like galaxies both the effects of general relativity and the role of a magnetic field can be ignored. However, in the 
case of active galactic nuclei with a huge central black hole and an important magnetic field, the Ernst spacetime can 
represent some relevant properties of the galactic structures. Therefore, it could even be astrophysically important to 
discuss and illustrate basic properties of the Ernst spacetime. 

The Ernst spacetime has the important property that it is not asymptotically flat. Far from the black hole, the 
spacetime is closely related to Melvin's magnetic universe, representing a cylindrically symmetric spacetime filled with 
an uniform magnetic field only. The Ernst spacetime is axially symmetric, and its structure corresponds to the structure 
of the Schwarzschild spacetime only along its axis of symmetry. Off the axis, the differences are very spectacular, and 
we shall demonstrate them for the equatorial plane, which is the symmetry plane of the spacetime. 

The line element of the Ernst spacetime reads 

ds 2 = A 2 |^-^l- — ^df 2 +^l- — ^ dr 2 + r 2 d0 2 

where M = M cgs G/c 2 is the mass, B = B^G 1 / 2 /c 2 is the strength of the magnetic field, A = 1 + B 2 r 2 sin 2 6. 

The dimensionless product BM 1 in astrophysically realistic situations. Really, BM ~ 2 x 10 _53 B C gs^cgs- 
For example, in case of a black hole with M — 10 9 M Q (typical for AGN), the value BM ~ 1 corresponds to 
B C gs ~ 10 11 Gauss, which is unrealistic. Some illustrative embedding diagrams are collected in Figure|4] 



r 2 sin 2 <9 
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■dcf> 2 , 



(22) 



FIGURE 4. Left column: ordinary geometry, right column: optical geometry of Ernst spacetime. It can be proved [ 5] that a critical 
magnetic field B c 0.0947 exists. For B > B c , neither throats nor bellies and no circular photon orbits exist. For B < B c , the throat 
and the belly develop, corresponding to the inner unstable and outer stable photon circular orbit. 

Kerr-Newman spacetimes 

We focus on the results concerning the optical reference geometry lT2lll . It can be shown [21] that the turning points 
of the embedding diagrams are really located at the radii where the centrifugal force vanishes and changes sign. Thus, 
we can conclude that it is exactly this property of optical geometry embeddings for the vacuum spherically symmetric 
spacetimes (see II HX |3l Q|]) survives in the Kerr-Newman spacetimes. However, photon circular orbits are displaced 
from the radii corresponding to the turning points of the embedding diagrams. 

After a straightforward, but rather tedious, discussion, the full classification can be made according to the properties 
of the embedding diagrams (see Clh . 

Examples of typical black-hole embedding diagrams of classes labeled in Stuchlfk et al. 02 111 as BH1-BH4 are 
presented in Figure [5] 




(a) (b) (c) (d) (e) 



FIGURE 5. (a), (b): Embedding diagram of the Kerr-Newman black holes of the type BHi , constructed for a 2 = 0. 16, e 2 = 0. 16. 
The rings in the 3D diagram represent photon circular orbits. Both corotating (gray spot in 2D diagram) and counterrotating (black 
spot in 2D diagram) orbits are displaced from the throat of the diagram where the centrifugal force vanishes. This is a general 
property of the rotating backgrounds, (c)-(e): Embedding diagrams of the black holes classes BH2-BH4. (c) Class BH2 (a 2 = 0.7, 
e 2 = 0.16); (d) class BH 3 (a 2 = 0.75, e 2 = 0.16); (e) class BH 4 (a 2 = 0.81, e 2 = 0.16). Taken from HQ]. 



CONCLUDING REMARKS 



Embedding diagrams of the optical geometry give an important tool of visualization and clarification of the dynamical 
behavior of test particles moving along equatorial circular orbits: we imagine that the motion is constrained to the 
surface z(p) (see yd). The shape of the surface zip) is directly related to the centrifugal acceleration. Within the 
upward sloping areas of the embedding diagram, the centrifugal acceleration points towards increasing values of r and 
the dynamics of test particles has an essentially Newtonian character. However, within the downward sloping areas 
of the embedding diagrams, the centrifugal acceleration has a radically non-Newtonian character as it points towards 
decreasing values of r. Such a kind of behavior appears where the diagrams have a throat or a belly. At the turning 
points of the diagram, the centrifugal acceleration vanishes and changes its sign. 
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